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I . INTRODUCTION 

For the determination of geopotential coefficients we can use data from ra- 
ther different sources, e.g. satellite tracking, gravimetry or altimetry. 

As each data type is particularly sensitive to certain wavelengths of the 
spherical harmonic coefficients it is of essential importance how they are 
treated in a combination solution. For example the longer wavelengths are 
well described by the coefficients of a model derived by satellite tracking, 
while other observation types such as gravity anomalies Ag and geoid heights 
N from altimetry contain only poor information for these long wavelengths. 
Therefore, the lower coefficients of the satellite model should be treated 
as being superior in the combination. In our contribution we present a new 
method which turns out to be highly suitable for this purpose due to its 
great flexibility combined with robustness. 

2 . METHODS 

In B . Middel/B. Schaffrin (1987) we introduced a method based on "robust 
collocation", which according to B . Schaffrin (1985,1986) is the Best homo- 
geneously Linear (weakly) Unbiased Prediction (hom-BLUP), as a promising 
technique for the combination of terrestrial gravity data with spherical har- 
monic coefficients from satellite tracking. With this method we obtain the 
predicted coefficients, collected in the vector x, by 

5 = (P g + P T ) _1 (P g ? s + P T S T a) (2.1) 

where contains the satellite coefficients and Pg is the corresponding 
weight matrix. Vector ^ contains coefficients which we obtained by a least- 
squares adjustment within a GauB-Markov Model from terrestrial gravity data 
and P^ is again the corresponding weight matrix. The terrestrial coefficient 
set is taken to be inferior with respect to the lower coefficents and by 
comparison with the satellite coefficient set £g we obtain the scalar factor 
a to fit to When the fitting factor is defined as a:=l we obtain the 

weighted mean oi both data sets being, of course, the "geodetic collocation" 
according to H. Moritz (1973) , which is the Best inhomogeneously Linear Pre- 
diction ( inhom-BLIP) . 

However hom-BLUP turned out to be robust against inconsistencies in and 
therefore, in this sense, superior to inhom-BLIP as we showed by applying 
statistical tests in B . Middel/B. Schaffrin (1988). Nevertheless we can 
make this approach more flexible by splitting up the weaker coefficient set 
into groups of a special character. We allow them specific fitting fac- 
tors collected in a vector a and name it "Mixed hom-BLUP", thus leading to 
the following solution: 
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( 2 . 2 ) 


- (Pg + P T ) ( p s ^ s + P T“t-^ ' 

In this formulation the vector £ is modified to a matrix ~ T where each col 
umn contains one group of coefficients at the respective places with zero 
entries otherwise. 


3. NUMERICAL EXAMPLE 


After this short description we now present some results of combination so- 
lutions up to degree 36 where we merged coefficients derived either from 

gravity anomalies or geoid heights or alternatively from a combination of 
them both, with GEM-L2 coefficients up to degree 20, as described m 
F J Lerch et a.1 (1982). For each combination we split up the weaker co 
efficient set £ by degree and by order when using the Mixed hom-BLUP tech- 
nique In figures 3.1 to 3.3 we plotted the components of the vector a 
with respect to the group (i.e. degree or order) together with their con- 
fidence intervals of a significance level of 95%. In addition, the scalar 
factors a of inhom-BLIP and hom-BLUP were added. 

Figures 3.1 show that we obtained with all methods different fitting factors 
a and therefore different results when we merge the GEM-L2 coefficients with 
coefficients £_ adjusted from gravity anomalies Ag. But it has to be men- 
tioned that whin using Mixed hom-BLUP and splitting up by order most of 
the components of a are not significantly different from the scalar a o 
ordinary hom-BLUP since these values lie inside the 95% confidence interva • 

The situation changes when we introduce coefficients £ T computed with geoid 
heights N. This data set is fully compatible with the GEM-L2 coefficien s 
and therefore all fitting factors are very close to 1 as illustrated in 
figures 3.2. In this case we obtained nearly identical solutions with al 

the methods. 

In figures 3.3 and 3.4 we present results which we obtained by a combination 
of GEM-L2 coefficients with coefficients which we adjusted from both data 
sets Ag and N. Figures 3.3 show that the scalar a is very close to 1 but 
with the Mixed hom-BLUP we obtained, in both cases of splitting, a-compo- 
nents far away from this value. Therefore, the solution with ordinary hom- 
BLUP leads to different results in both cases of splitting up the coeffi- 
cients £ 


This can be very clearly seen in figure 3.4 where the degree 
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variances of the different solutions have been plotted. 


4. 


CONCLUSIONS 


The use of "Mixed hom-BLUP", which we present in this cont |; l ^ ut ’; on ’ 
to different solutions compared with "geodetic collocation (mhom BL ) 
and "robust collocation" (hom-BLUP) if the combined data sets are not com- 
patible. Due to the great flexibility and robustness of this metho we 
expect that it is highly suitable for estimating geopotential coefficients 
when combining heterogeneous data sets. 
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